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Abstract. In 5 Scuti star models, the calculated amplitude ratios and phase differences for multi-colour 
photometry exhibit a strong dependence on convection. These observables are tools for determination of the 
spherical harmonic degree, I, of the excited modes. The dependence on convection enters through the complex 
parameter /, which describes bolometric flux perturbation. We present a method of simultaneous determination 
of / and harmonic degree t from multi-colour data and apply it to three 5 Scuti stars. The method indeed works. 
Determination of £ appears unique and the inferred f's are sufficiently accurate to yield a useful constraint 
on models of stellar convection. Furthermore, the method helps to refine stellar parameters, especially if the 
identified mode is radial. 
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1. Introduction 

The 8 Scuti stars are pulsating variables located in the 
HR diagram at the intersection of the classical instability 
strip with the main sequence and somewhat above it. The 
observables of primary interest for asteroseismology are 
oscillation frequencies. However, information about am- 
plitudes and phases of oscillations in various photometric 
passbands is also useful. So far the main application of 
multi-colour photometry of 8 Scuti stars has been deter- 
mination of £ degree of observed modes (see e.g. Balona 
& Evers 1999, Garrido 2000). This is an important ap- 
plication because knowledge of £ is an essential step for 
mode identification. The £ diagnostic makes use of dia- 
grams in which the amplitude ratio determined in two 
passbands is plotted against the corresponding phase dif- 
ference. The observational data are compared with ranges 
calculated for relevant stellar models and assumed I val- 
ues. The ranges reflect uncertainties in stellar parameters 
and physics. The inference is easy if the ranges do not over- 
lap. This is largely true for (3 Cephei stars, but not for 8 
Scuti stars. For the latter objects, a major uncertainty in 
the calculated ranges arises from lack of adequate theory 
of stellar convection. 
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Calculations of the amplitude ratios and phase differ- 
ences make use of the complex parameter /, which gives 
the ratio of the radiative flux perturbation to the radial 
displacement at photosphere. The parameter is obtained 
with the linear nonadiabatic calculations of stellar oscil- 
lations. The problem, which has been already emphasized 
by Balona & Evers (1999), is that / is very sensitive to 
convection whose treatment still remains rather uncertain. 

The strong sensitivity of calculated mode positions in 
the diagnostic diagrams to treatment of convection is not 
necessarily a bad news. Having data from more than two 
passbands we may try to determine simultaneously £ and 
/. If we succeed, the /-value inferred from the data would 
yield then a valuable constraint on models of stellar con- 
vection. The aim of our work is to examine prospect for ex- 
tracting / from multi-colour photometry of 8 Scuti stars. 

In the next section we demonstrate strong sensitivity 
of / and, in a consequence, of mode position in the diag- 
nostic diagrams, to the description of the convective flux. 
Our treatment of convection is very simplistic. We rely 
on the mixing-length theory and the convective flux freez- 
ing approximation. We study how calculated positions in 
the diagnostic diagrams vary with the changes of mixing- 
length parameter, a. 

Our method of inferring / from the data is described 
in Sect. 3. The subsequent section presents application of 
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Fig. 1. The real and imaginary part of the /-parameter for radial oscillation of a 1.9 M & star in the main-sequence 
phase. In the panels from left to right we show the /-parameter for three values of a. Dots represent the unstable 
modes, whereas open circles the stable ones. 



the method to several S Scuti stars. In the last section we 
summarize results of our analysis. 

2. Calculated mode positions in the diagnostic 
diagrams 

We use here the standard description of oscillating stellar 
photospheres (cf. Cugier et al. 1994). The local displace- 
ment is adopted in the form 

Sr(R, 6, if) = eRRe{Y e m e- iuJt }, 

where e is a small complex parameter fixing mode ampli- 
tude and phase. The associated perturbation of the bolo- 
metric flux, Thoi, and the local gravity, g, are then given 

by 



= eRe{/^ m e- lwt }, 



and 



^ = -12 
9 



co 2 R 3 \ SR 
GM ) R' 



With the static plane-parallel approximation for the 
atmosphere we can express the complex amplitude of 



the monochromatic flux variation as follows (see e.g. 
Daszyriska-Daszkiewicz et al. 2002) 



A\i) = eYT(i, 0)b$(Dt t + D 2 . e + D* e ), 



(1) 



where i is the inclination angle and A identifies the pass- 
band, 



1 d\Qg{Ty\b\\) 

A J aiogr eff : 



D 2 , t = (2 + £)(!-£), 



D 



3,e 



uj 2 R 3 



dlog(T x \b}\) 



GM ' J dlogg ■ 
The disc averaging factor, bi, is defined by the integral 



b$ = 



h\(/j,)(j,Pe{[i)d(i, 



where the function h\ describes the limb darkening law 
and Pi is the Legendre polynomial. The partial deriva- 
tives of J-\\bg I may be calculated numerically from tabular 
data. Here we rely on Kurucz (1998) models and Claret 
(2000) computations of limb darkening coefficients. With 
Eq. 1 we can directly obtain the amplitude ratio and the 
phase difference for chosen pair of passbands, which are 
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called nonadiabatic observables. Nonadiabaticity of oscil- 
lations enters through the complex parameter /, which is 
the central quantity of our paper. 

In Fig. 1 we illustrate how the choice of the mixing- 
length parameter, a, affects the value of /. We can see 
the large effect of the choice, particularly between a = 0.0 
and a = 1.0 in the cooler part of the sequence where all 
modes are unstable. We have found that important is only 
the value of a in the H ionization zone. Models calculated 
with a fixed in this zone and varied in the Hell ionization 
zone yield very similar values of /. 

In the next two figures we show how the differences 
in a are reflected in the Aj,- V /A y vs. ifib-y — fy diagram 
employing the Stromgren passbands . The effect of varying 
a is very large indeed. First, in Fig. 2 we show all the 
unstable modes in the frequency range covering pi to p 4 
radial mode. The domains of different £'s partially overlap 
even at specified a. The ambiguity is removed once model 
parameters are fixed, as seen in Fig. 3, but sensitivity to 
a is well visible. The sensitivity indicates that, if we are 
able to deduce values of / from multi-color photometry, 
we will have at hand a valuable new constraints on stellar 
convection. 

The test should be applied to more realistic modeling 
of convection - pulsation interaction than used here. Our 
only aim here was to show sensitivity of the nonadiabatic 
observables to convection and we believe that our approx- 
imation is adequate for this aim. In section 4 we compare 
the calculated and measured /'s. 

3. A method for inferring / values from 
observations 

We begin with rewriting Eq. 1 in the form of the following 
linear equation. 



(2) 



where 



e = eY e m (i,0) 
= b x (D 2 



h x n x 



Here the A superscript identifies the passband. Eqs.(2) for 
a number of A's form a set of observational equations. 
In the right-hand side we have measured amplitudes, A x 1 
expressed in the complex form. The quantities to be de- 
termined are (if) and e. Both must be regarded complex. 
Of course we are primarily interested in the value of /. 
However, inferred value of e may also be useful as a con- 
straint on mode identification if it is found unacceptably 
large. 

Having data only for two passbands we can infer / in 
a unique way once we know t. However, we usually do not 
and therefore we need at least three passband data. The 
procedure is to determine / by means of x 2 minimization, 
assuming trial values of £. We will regard the I and the 
associated complex / value as the solution if it corresponds 
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Fig. 2. The effect of mixing-length parameter on the lo- 
cations of modes with different £ degree. Here we show 
the positions of £ = 0, 1, 2 unstable modes in the diagnos- 
tic diagrams involving b — y and y Stromgren filters for 5 
Scuti models of 1.9M©. The frequency range covers pi to 
Pi radial modes. The upper and lower panels are for a = 
and a = 1.0, respectively. 



to x 2 minimum which is significantly deeper than at other 
■f s. We will consider the £ values up to six. 

If we have data on spectral line variations, the set of 
equations (2) may be supplemented with an expression 
relating e to complex amplitudes of the first moments, 
M x , 

XjjR[u ) + 9M±\ i=M x . (4 ) 



where 



v\ = £\ h x {\i)\i(Pi-x(ii) - fiPt(ji))dn. 
Jo 

are coefficients representing a straightforward general- 
ization of the coefficients ug and vg introduced by 
Dzicmbowski (1977) for gray atmospheres. We stress that 
the first moment is the only measure of the radial velocity 
amplitude, which does not depend on the aspect and the 
azimuthal number, m, just like the light amplitude. 
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Table 1. Stellar parameters for three 6 Set stars: (3 Cas, 20 CVn and AB Cas. 



Object 


Sp 


Period [d] 


7r [mas] 


logT off 


logL 


logs* 


log g * * 


[m/H] 


v e sin i 


P Cas 


F2III-IV 


0.1009 


59.89 ± 0.56 


3.856 ± 0.01 


1.426 ± 0.008 


3.60 ±0.01 


3.66 


0.0 


69 


20 CVn 


F3III 


0.1217 


11.39 ±0.69 


3.874 ±0.01 


1.881 ± 0.053 


3.62 ±0.01 


3.49 


0.5 


15 


AB Cas 


A3V+KV 


0.0583 


3.33 ±1.30 


3.908 ±0.01 


0.998 ± 0.430 


4.55 ±0.01 


4.26 


0.0 


55 
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Fig. 3. The same as in Fig. 2, but the one stellar model 
with logoff = 3.867 is considered. The upper panel shows 
the photometric observables calculated at a = 0.0, and the 
lower one those at a = 1.0. 



There are uncertainties in model parameters, which 
enter the expressions for D x and . The partial deriva- 
tives of .Fa 1 6^ |, which appear there depend on T c g, logy 
and the metallicity parameter [m/H]. We will not con- 
sider here the stars with chemical peculiarities, thus we 
are adopting the solar mixtures of heavy elements. We 
repeat our minimization for the ranges of these three pa- 
rameters consistent with observational errors as well as 
with our evolutionary tracks. 

If the best fit is for I = 0, an additional constraint on 
models follows from mode frequency, because the I = 
frequency spectrum is sparse and the radial order of the 
mode may be easily identified. In this case we can tune a 
star in the exact values of the observed frequency. Stellar 
parameters have significant effect on \ 2 - We will see that 



based on we can obtain more stringent constrains on 
these parameters. 

Generalization of the method to the case of modes cou- 
pled by rotation is straightforward. Such modes are repre- 
sented in terms of a superposition of spherical harmonics 
with ^'s differing by 2 and the same m's. The monochro- 
matic amplitude of a coupled mode is given by 



where a>k coefficients are solutions of the degenerate per- 
turbation theory for slowly rotating stars (see Daszyhska- 
Daszkiewicz et al. 2002). The values of a k coefficients de- 
pend on the rotation rate and the frequency distance be- 
tween modes. 

The counterpart of Eq.(2) is 



^w x k [Vl(ef) 



k 



■£le] = A> 



(5) 



where 



w X k =Y%{i,Q)bU k . 



The main difference relative to the single £ case is that 
now the result depends on the aspect. Further, we expect 
a strong dependence of calculated amplitude on model 
parameters because the values of a k are very sensitive to 
small frequency distances between coupled modes. 

4. Applications 

We applied the method described above to data on three 
5 Set variables: (3 Cas, 20 CVn and AB Cas. In Table 1 
we give parameters for these stars. In this table we rely 
mostly on the catalogue of Rodriguez et al. (2000). 

The photometric data were dereddened according to 
Crawford (1979) and Crawford & Mandwewala (1976), 
and then the effective temperatures, gravities and bolo- 
metric corrections were obtained from Kurucz's (1998) 
tabular data. In Table 1 we give two values of logg, from 
Kurucz data and from evolutionary tracks. To calculate 
the D x and coefficients we took the second one, which 
we regard more reliable. Errors in logT e ff and log g corre- 
spond to typical errors from the photometric calibration 
procedure. 

For (3 Cas and AB Cas we adopted solar metal abun- 
dance, whereas for 20 CVn we used [m/H] = 0.5 (Z ps 
0.06) due to Hauck et al.(1985) and Rodriguez et al.(1998) 
and checked also the other one, [m/H] = 0.3 (Z w 0.04). 
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Fig. 4. In the left panel we show the x 2 as a function of I obtained for three models of (3 Cas star. In the right panel 
we plot the x 2 as a function of inclination considering the mode coupling between I — and I — 2 at v YO t =70 km/s. 



The values of log L for (3 Cas and 20 CVn are derived from 
the Hipparcos parallaxes. The parallax for AB Cas is very 
inaccurate. The central value locates this object close to 
the ZAMS. Though the star is a component of an Algol- 
type binary system, it is reasonable to assume that it is 
described by ordinary mass-conserving models because it 
was originally less massive and the episode of the rapid 
mass accretion is most likely forgotten. This is what we 
assume in this paper. Therefore, we adopted as the min- 
imum luminosity the value at ZAMS, and the maximum 
luminosity as the highest values allowed by the Hipparcos 
parallax. 

4.1. (3 Cas 

(3 Cas is a S Set star with radial velocity variations of 
2 km/s (Mellor 1917). Mills (1966) classified it as S Set 
variable on the basis of photometric observations. 

f3 Cas is one of a few S Set stars in which only one 
mode has been detected so far. Rodriguez et al. (1992) 
identified this mode as P2 or p^ of £ — 1, on the basis 
of photometric diagrams in Stromgren filters. Balona & 
Evers (1999), using the same photometric data, did not 
get an unambiguous identification. 

We have rather precise parameters for this, the bright- 
est and the nearest S Scuti variable. Note, in partic- 
ular, small luminosity errors in Tab. 1. The adopted 
metal abundance, Z = 0.02, is based on the [m/H] 
value obtained from IUE spectra by Daszyhska & Cugier 
(2002). The photometric amplitudes and phases are from 
Rodriguez et al. (1992). 

The frequency value combined with mean density im- 
plies that if the mode were radial, it could be only p^. 
We first assume that the mode is adequately described in 
terms of single spherical harmonic and consider ^-values 
from to 6. The uncertainties in the values of M and R 



10 



6 



I 


1 I 1 I 


1 i 1 i 1 i 1 i 1 


i > i 














q \^a=0.5 








c 










a=1.0 








s 






o log r eS = 3.856 








* log T =3.846 f 








cx s4a=1.6 


I 


i . i 


a=2.5 




i . i 



-8 -6 -4 -2 2 4 6 8 



Fig. 5. Comparison of the / values inferred from 
Stromgren photometry for (3 Cas with the values from the 
models calculated with various values of the MLT param- 
eter a and the mass of 1.95-Mq. Points with the error bars 
represent inferred values assuming two indicated effective 
temperatures. 



are inconsequential and we consider only uncertainty in 
T e ff. The estimated value of the mass is 1.95M©. In the 
left panel of Fig. 4 we see that the uncertainty on the ef- 
fective temperature does not impair identification of the 

1 value. The x 2 minimum at I = 1 is the deepest one, 
particularly at the two lower effective temperatures. 

The star is a relatively rapid rotator, w rot is at least 70 
km/s, therefore we have to consider the possibility that 
the mode is a coupled one, most likely an t = and £ = 

2 superposition (see Daszynska-Daszkiewicz at al. 2002). 
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Fig. 6. Dependence of x 2 on (- f° r f° ur models of 20 CVn 
on the edges of the error box obtained for Z = 0.06. 

Relying on the formalism described at the end of Section 
3, we evaluated the best /-values and assiociated \ 2 as 
a function of the aspect. In the right panel we show \ 2 
variation with the inclination. We can see that at no value 
of the inclination the \ 2 is nearly as low as at £ = 1 . Thus 
we conclude that the mode excited in (3 Cas is most likely 
a single £ — 1 mode. 

Values of / corresponding to this identification are 
shown in Fig. 5 together with the error bars, which are 
the errors from the least-square method. The uncertainty 
in temperature is more significant than the errors of / de- 
termination. The theoretical / values were calculated for 
M — 1.95M© assuming five values of the mixing-length 
parameter a: 0.0, 0.5, 1.0, 1.6 and 2.5. Still the constraint 
on convection is interesting because the range of the ac- 
ceptable values of / is narrower than the range of the 
calculated values with different a's. 

The observed values of fn are closer to those calcu- 
lated with a = 0, which may be taken as an evidence that 
convection in the H ionization zone is relatively inefficient. 
However, values of // require rather higher a's (about 1) 
instead. In view of our crude treatment of convection- 
pulsation interaction we have to take these indications 
with a great caution. 

4.2. 20 CVn 

This S Scuti variable is also regarded to be monoperiodic 
(e.g. Shaw 1976, Pena & Gonzalez-Bedolla 1981). The 
mode was identifited as £ = by means of photometry 
(Rodriguez et al. 1998) as well as spectroscopy (Chadid et 
al. 2001). Our result shown in Fig. 6 for Z = 0.06 clearly 
confirms the previous identification. For models on the 
edges of the error box the minimum of x 2 is at £ — 0. The 
same is true for Z = 0.04, but the values of the \ 2 are 
higher. 

As we mentioned in Sect. 2, having such identification 
of pulsating mode we can refine stellar parameters by fit- 




3.92 3.90 3.88 3.86 3.84 3.82 

lo 9 r 8 „ 

Fig. 7. The observational error box for 20 CVn on HR 
diagram. The lines of constant radial order, n = 3, 4, 5, 
are also drawn. 




Fig. 8. The variation of \ 2 with effective temperature for 
the models along the lines n = 3 and n = 4. 



ting the observed period. Still, we have to consider various 
radial orders, n. In Fig. 7 we show the HR diagram with 
the error box representing uncertainty of stellar parame- 
ters and the lines of the constant period (P = 0.1217 d) for 
n = 3, 4, 5, obtained with evolutionary models calculated 
for Z = 0.06 and indicated masses. Only models along 
these lines are allowed. We can see that as for the radial 
order we have only two possibilities, n — 3 or n — 4. In 
the case of Z = 0.04 only n — 4 is allowed but the \ 2 is 
significantly larger. We thus see that our method allows 
to constrain the Z value. 

In Fig. 8 we show the variations of the x 2 with logT c g 
for the tuned models with n = 3 and n = 4. The vertical 
lines correspond to the intersections of the constant period 
lines with the error box, so that only the values of log T c s 
between these lines for a given n are allowed. The models 
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Fig. 9. The same as in Fig. 5, but with empirical values 
of / for the models with minx 2 f° r U = 3 and n = 4 of 20 
CVn. 

yielding the lowest x 2 are just models 1 and 3 in Fig. 7 for 
n = 4 and n = 3, respectively. The values of % 2 are 0.19 
for n = 3 and 0.24 for n = 4, what indicates that both are 
possible as identifications of the radial order in this star. 

Fig. 9 shows the empirical values of the nonadiabatic 
/-parameter for models corresponding to the deepest x 2 
minima for four values of a: 0.0, 1.0, 1.6 and 2.5. Relative 
positions of empirical and theoretical / values in 20 CVn 
are qualitatively similar to those in (3 Cas (cf. Fig. 5). 
However, these values, both calculated and inferred, are 
considerably higher than in the case of /? Cas, which is a 
consequence of higher radial order. The I — 1 mode iden- 
tified in (3 Cas has frequency between the n — 2 and n = 3 
radial modes. 

4.3. A B Cas 

The next example is the primary component of an Algol- 
type system. As we have already pointed out, luminosity 
of this star is very uncertain due to the large error in the 
parallax. We will see that with our method we can signif- 
icantly improve the accuracy of the stellar parameters as 
provided that the identified mode is radial. 

In the whole range of allowed parameters % 2 has by far 
the deepest minimum at I = 0, as we can see in Fig. 10. 
This identification is in agreement with that of Rodriguez 
et al. (1998). In Fig. 11 we show the HR diagram with 
the error box representing uncertainty of stellar param- 
eters and the lines of the constant period (P = 0.0583 
d) for n — 1,2,3, obtained with evolutionary models cal- 
culated for Z = 0.02 and indicated masses. The deep- 
est minimum of x 2 (0.25) along the n = 1 line occurs 
at M = 1.77, log Teg = 3.8985, logL = 1.047. The lowest 
value of x 2 for n = 2 (0.62) is at M = 1.91, logT cff = 
3.8985, logL = 1.217, and for n = 3we got X 2 = 1-02, 



x 




I I i I i I i I i I i I i L 

1 2 3 4 5 6 



I 

Fig. 10. Dependence of % 2 on I for four models on the 
edges of the error box of AB Cas. 
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Fig. 11. The observational error box for AB Cas on HR 
diagram. The lines of constant radial order with n — 1, 2, 3 
are also shown. 

T e g = 3.8985, logL = 1.361. We can see that n = 1 iden- 
tification is strongly favored. For such mode identification 
the value of logL is between 1.047 and 1.144 leading to 
X 2 = 0.25 and 0.56, respectively. 

In Fig. 12 we compare empirical values of / with ones 
calculated for various mixing-length parameter, a. The 
star is the hottest of the three objects. Still, calculated 
values are strongly affected by the choice of a. Just like 
the previous, regardless the value of a, the empirical val- 
ues of the real and imaginary parts of /, could not be 
simultaneously reproduced with our model calculations. 

5. Conclusions and discussion 

We believe that the three examples of 8 Scuti stars consid- 
ered in the previous section, clearly show that there is a 
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Fig. 12. The same as in Fig. 5, but with empirical values 
of / for two tuned models of AB Cas with n = 1 and 
logT eff = 3.908 and 3.898. 

wider application of multicolor data on excited mode than 
explored so far. In addition to the spherical harmonic de- 
gree, £, we were able with our x 2 minimization method to 
infer the value of the complex parameter / which describes 
the bolometric flux perturbation and yields a strong con- 
straint on models of stellar convection. In the two cases 
when determined value of I was zero, we were able to de- 
termine the radial order of the mode and significantly limit 
the uncertainty of the stellar parameters. 

The values of / we found could not be reproduced with 
model calculations made with the convective-flux-freezing 
approximation for any value of the MLT parameter a. 
We should not be surprised. After all, the approximation 
adopted is grossly inadequate. Though there are certain 
common features in the relative positions of the deduced 
and calculated /'s in the three cases considered, it is pre- 
mature to draw any conclusions about properties of stellar 
convection from this fact. 

Inadequacy in our treatment of convection is not the 
only possible cause of the large discrepancy between the 
calculated and empirical /'s. Our use of the Eddington ap- 
proximation in calculation of the nonadiabatic oscillations 
may be an oversimplification. We believe it is of secondary 
importance as the value of / is determined at the optical 
depth t 1. Our use of static atmospheric models with 
the depth-independent g seems also a good approxima- 
tion. Still, it should be kept that these two approximation 
should be at some point verified. The accuracy of the at- 
mospheric models is of a greater concern, which was dis- 
cussed recently by Heiter et al.(2002), where sensitivity of 
the atmospheric structure and observable quantities to the 
convection treatment was demonstrated (see also Smalley 
& Kupka 1997). Whatever is the cause of discrepancy, the 
goal in model calculations should be to achieve consistent 
f's. 



Determination of £ is an independent goal. It is impor- 
tant that it could be done without an a priori knowledge 
of /. This is not a new finding. This has been done earlier 
for both 5 Scuti and f3 Cephei variables. In particular, in 
the latter case, when / may be approximately treated as 
real, the three color data often allow for an unambiguous 
t determination (e.g. Heynderickx et al. 1994, Balonafe 
Evers 1999 ). For 5 Scuti stars the situation is more com- 
plicated. Although, we believe, our method represents an 
improvement, it still does not always work . We may be- 
lieve in the inferred values of both / and I only if the 
minima of \ 2 are strongly I dependent. This was the true 
in all the three cases considered in the previous section. It 
is not always so. For the 5 Scuti star 1 Mon, for instance, 
we found \ 2 quite flat between £ = and 2. For this star 
we also made an unsuccessful attempt to use the data 
of the radial velocity as provided by Balona et al.(2001), 
but even including this quantity did not change the re- 
sult. This failure should not discourage future efforts to 
combine spectroscopic data. 

The main message of our work is that applications of 
multi-color photometry data of 5 Scuti stars go beyond 
identification of the £ values of the excited mode. We 
showed that such data allow to refine parameters of the 
stars and, what we regard most important, yield strong 
constraints on models of stellar convection. 
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